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LONG-TERM  GOALS 

The  long  tenn  goal  of  this  project  is  to  advance  the  state  of  the  art  in  mathematical  methods  for 
detecting  key  transport  structures  in  velocity  field  data  sets  for  spatially  complex,  time-dependent, 
ocean  surface  flows.  Such  transport  structures  are  typically  not  inherently  obvious  in  snapshots  of  the 
Eulerian  velocity  field  and  require  analysis  methods  with  foundations  in  Lagrangian  flow  transport. 
Within  the  context  of  this  DRI,  a  long  tenn  goal  is  to  improve  understanding  of  the  formation  and 
location  of  submesoscale  fronts  in  the  South  China  Sea  (SCS). 

OBJECTIVES 

The  scientific  objective  is  to  test  and  develop  novel  methods,  with  a  focus  of  the  approach  of 
Lagrangian  Coherent  Structures  (LCSs).  We  seek  to  demonstrate  the  utility  and  robustness  of  the 
method  through  application  to  ocean  data  sets  pertaining  the  SCS,  and  more  specifically  the  region 
immediately  to  the  southwest  of  Taiwan.  Furthermore,  we  seek  to  directly  apply  these  methods  to  data 
sets  from  HF  radar. 

APPROACH 

The  approach  involves  using  MATLAB  and  C++  based  software  to  perform  processing  of  the  data 
sets.  The  code  has  been  developed  by  a  graduate  student  (Allshouse)  and  a  postdoc  (Leclair).  Data  sets 
for  analysis  are  provided  by  other  members  of  the  DRI,  such  as  John  Wilkin,  and  HF  radar  data  from 
the  Taiwanese  TORI  network. 

WORK  COMPLETED 

There  have  been  substantial  developments  of  the  computational  code.  First,  and  foremost,  the  code  has 
been  reconfigured  to  work  in  arbitary  curvilinear  coordinates.  This  has  several  benefits.  For 
investigations  of  large  scale  ocean  domains  the  LCS  code  can  work  directly  with  the  girdded 
information  provided  by  a  numerical  model  and  no  longer  approximated  the  domain  by  projection  on  a 
rectilinear  grid.  This  also  reduces  errors  introduced  into  the  trajectories  of  advected  material.  Finally, 
in  these  native  coordinates  the  LCS  code  is  more  capable  of  handling  the  complex  shapes  of  coastlines, 
which  is  an  issue  that  is  known  to  introduce  suprious  LCS.  The  second  advance  of  the  code  is  that  a 
scheme  has  been  introduce  to  ensure  conservation  of  flux  through  a  grid  cell,  so  that  there  is  no 
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numerical  introduction  or  extraction  of  volume,  again  creating  spurious  LCS.  Finally,  the  code  has 
been  configured  to  run  on  a  GPU  so  that  it  runs  considerably  faster. 

In  February  2014,  the  postdoc  working  on  this  project  participated  in  the  South  China  Sea  field 
experiments  and  also  visited  TORI.  The  HF  radar  data  from  TORI  was  processed  in  near  real  time  and 
supplied  to  the  science  team  in  order  to  facilitate  their  search  for  submesoscale  fronts  in  the  region. 

RESULTS 

The  LCS  processing  during  the  SCS  cruise  was  able  to  detect  Lagrangian  features  characteristic  of 
fronts  and  relay  this  information  to  cruise.  An  example  of  this  processing  applied  to  a  TORI  data  set  is 
presented  in  figure  1 .  There  were  challenges  in  the  quality  of  the  data  due  to  the  time-varying  technical 
issues  with  the  domain  of  coverage  of  the  HF  radar.  More  work  needs  to  be  done  to  detennine  if  the 
LCS  detected  in  the  HF  radar  data  truly  corresponded  to  ocean  fronts  detected  by  the  science  team. 
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Figure  1:  Attracting  LCS  for  a  12  hour  time  window  in  February  2014  using  TORI  HF  radar 
data.  Several  attracting  features  are  observable  and  are  candidates  for  surface  signatures  of 

ocean  fronts. 

We  have  also  investigated  the  application  of  the  methods  to  the  study  of  regions  containing  coral  reefs, 
which  could  be  applicable  to  the  vicinity  of  Dongsha  Island,  for  example.  Figure  2  shows  the 
connectivity  of  a  coral  reef  system  for  a  12  hour  time  window  calculated  using  the  LCS  approach. 
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Figure  2:  (top  left  and  right)  Repelling  and  attracting  LCS  for  a  coral  reef  system  (off  the  west  coast 
of  Australia)  over  a  12  hour  time  window,  (bottom  left  and  right)  The  flow  domain  is  organized  into 
three  components:  offshore  vortices  that  trap  fluid  (yellow),  northward  advected  material  (magenta) 

and  southward  advected  material  (green). 


IMPACT/APPLICATIONS 

The  analysis  during  the  SCS  experiment  provided  valuable  guidance  for  the  real-time  decision  making 
of  the  science  team.  The  improved  computational  methods  make  the  LCS  approach  ever  more 
computationally  efficient  and  robust,  further  advocating  its  real  world  applicability.  A  major  LCS 
workshop  was  held  at  the  Banff  International  Research  station. 


3 


